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We investigate the QCD phase diagram as a function of isospin chemical potential at a fixed 
temperature by directly putting large numbers of tt+s into the system. Correlation functions of 
7r+-systems involves A^!A^! contractions, and become extremely expensive when is large. 
In order to alleviate this problem, a recursion relation of correlation functions has been derived 
in Ref. [1] that substantially reduces the number of independent contractions needed and makes 
the study of many pions systems be possible. In this proceeding this method is investigated 
numerically. We have also constructed a new method that is even more efficient, enabling us to 
study systems of up to 72 tt+'s. 
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1. Introduction 

The QCD phase diagram has been studied with many methods, but there are still many ques- 
tions waiting to be answered. At zero isospin density, there is a crossover (phase transition) from 
a confinement phase at low temperature to a decomfinement phase at higher temperature. From 
chiral perturbation theory (;t;PT), we expect that when the isospin chemical potential reaches the 
mass of a single pion, pions start to condense, and at higher isospin chemical potential there might 
be a crossover or a phase transition to another state [2]. Systems containing up to 12 tt+s have been 
studied by directly computing the required contractions in Ref. [5] [3], finding agreement with the 
expectations from order to put more TT+s into the system, a second source is required be- 

cause of the PauU principle, and studying systems of more than 12 71+ becomes almost impossible 
by directly computing all possible contractions. 

In Ref. [1], a recursion relation of correlation functions of «-7r+ systems is constructed which 
significantly reduces the amount of work needed, enabling us to study systems up to 24 ;r+'s. 
Recursion relations of M species systems have also been discussed in the same paper. Even with 
the recursion relation, studying a 3-source system (36 pions) still consumes a substantial amount of 
time. In order to overcome this problem we have constructed a new method which allows practical 
calculation of up to 727r+s. 

In this proceeding, we first review the methodology of the recursion relation in position space 
and then extend its application to momentum space, after that we present studies of systems con- 
taining up to 247r+s by applying the recursion relation in momentum space. At last we present 
results from the new method to study large number of pion systems. 



2. Methodology of the recursion relation 

2.1 Recursion relation in position space 

In order to explore systems containing up to 12M 71+ 's, M different sources are required be- 
cause of the Pauh principle. A correlation function of a system with rij-n^s in the /''^ source is: 

C(„,<,...,„„,.,^)(0 = ^ (l ?r+(x,o) (^-{yi,0)^ '...(7r-(y„„0)) "' ^ , (2.1) 

where fi = Calculating this correlation function according to Wick's theorem involves 

h\h\ contractions, which make the study for a system of large number of 7r+'s extremely time 
consuming. However the recursion relation of uncontracted correlation functions 2(„i,„2,...,,j,„)(0' 
makes the study of such systems feasible. The correlation function can be expanded as: 

n,„n-)it) = n ) ( Q(nun....,n„,){t) > , (2.2) 

where( ) denotes spin color trace, and 2(„i,„2.. ..,«,„) (0 satisfies the following ascending recursion 
relation: 

Q{ni + l,n2,...,n,„) = ( 2(«i,;i2, ...,«„)) Pi — n Q[nx,n2,...,n„,)P\ 

••• + ( 2(ni + l,n2,---n*-l, ■-,"»,)) Pk ~ ^ Q{ni + \,n2,...nt-\,...,nm)^k 

■■■ + {Q(rn + \,n2,...,n„,-\)) Pm — ^ Q(ni + \.n2,...n„-\)Pm , (2.3) 
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and the initial conditions are 2{i,o,...,o) = A = > 2(o,i,...,o) = Pi = ^2, • " " > where A;'s are uncon- 
tracted single pion correlators. Descending recursion relations can also be derived: 



M 

e„ = I 
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(2.4) 



where n = {ni,n2,--- ,«m)> 1/t = (0,0,--- ,1,0,---) with only the k'^ nonvanishing unit element, 
and 2i2,...,i2, Pk and A are constructed in the following way: 

ei2,...,i2 = {N-l)\det{A)-lN 
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(2.5) 



where Aij{t) is defined as Ajj (t) = Y,xS (x;, x) 5+ (xy, x). 

One way to construct A, j is shown in the left plot of Fig. 1 . Correlation functions of two 
species from multiple sources have similar recursion relations, for detail discussion about the re- 
cursion relations see Ref. [1]. 

2.2 The recursion relation in momentum space 

A correlation function of a system having ?ii-7r+s in the first source and n2-7l^s in another 
source with total momentum n\Pf^ +n2P/2 is: 

C„.;,.,„,.. it) = (tl(^L ^'-'(P'''"-P^'"')^^(x,-,075^(x/,-,0^ • n (l,e'^'j''d{yj,0)r5uiyj,0)^ ^ 

where n = n\ + n2. Momentum conservation requires that nip} + «2Pi —n\p\ — niPi = Ly=i P/; 
must be satisfied to get non- vanishing C„^jc+^„^j^+. Each choice of P'j,i,j = li2 satisfying this 
relation is a separate measurement. By replacing propagators in position space by propagators in 
momentum space, a similar recursion relation still holds. The only difference is the construction 
of uncontracted correlation functions Aij defined as A,j- (t) = Lp*^ (P/ )P) (pj^P — Pf)^^ which 
are compared on Fig. 1. 



3. Results 

Because of the finiteness of the temporal extent and the factorisable nature of the multi-hadron 
systems being constructed, thermal effects are particulai^ly important in multi-hadron systems and 
the correlation functions assume the form: 

Cn^-it) = t ('')z:e-^''"'+^"-"'^^^^cosh{{E,n-E„^^)-{t-T/2)) + --- (3.1) 

m=0 \ / 

where dots represent higher excitations, T is the temporal extent and is the energy of a rest 
system of «-7r+s. The ground state comes from all 71+ 's propagating in the same direction in time, 
and thermal states are from some 71+ 's propagating in one direction while the rest propagate in the 
opposite direction. 
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(a) (b) 



Figure 1: The left figure shows how to construct uncontracted correlators in spatial space and the right figure 
shows how to get its counterparts in momentum space. A, ^ is constructed by following the line from source 
/ to x{p), returning to source /', multiplying the S{S' ) with respect to each line and summing over x{p). 



3.1 Verify the dispersion relationship 

Calculations have been performed using anisotropic { 16^ , 20^ , 24^} x 128 lattices with anisotropy 
^ = 3.5 at a quark mass corresponding to a m^j = 390Mev. On the lattice, only discrete momentum 

are allowed. E,„-^+ of systems with total momentum p, = m-p, forp = (0,0, 1), (0, 1, 1), (0,0,2), 

have been extracted for m = 4, 3, 2 respectively and are fitted into the dispersion relation: — 2^ ~ 

(^)^ ~ ^ m"'**^ > where p, = m-p returning |c| = 1.015(32), which confirms that these many 
hadron systems are describable in special relativity (thermal states would have a different dependence). 

3.2 One species from single source 

For notational convenience, in the following p} is denoted as pi , P2 as p2, Pj as p3 and as p4. 
Azimuthal symmetry ensures many combinations of pi, p2 be separate measurements of the same 
physics, which provide more statistics. As £'„;j+ extracted from different choices of pi and p2 agree 
with each other within errors, we choose pi = p2 = ( 1 , 1 , 1 ) for further discussions. Extracting E„j,^+ 
by fitting to functions with and without one excited state in addition to the thermal states discussed 
above from different time intervals produces consistent results, as is shown in Fig. 2. 

Decomposing C„ji+{t) into different contributions gives insight into how much each state in 
Eq. 3.1 contributes. In Fig. 2, C\2n+{t) is shown. In this figure, the green line is from the first 
excited state, the blue one is the ground state and other lines are thermal states. It is remarkable 
that the zero temperature ground state is not dominant in any region of this correlator. Even for 
Cnn+ (0 ground state is strongly contaminated by both thermal states and excited states, and it 
becomes difficult to extract the ground state energy from €„■„;+ (t), n > 12, with the temporal extent 
of these configurations. 

3.3 One species from two sources 

By choosing pi = p2, P3 = P4 but pi 7^ p3 coiTclation functions of systems having up to 24?!+ 
have been computed with the same recursion relation, but E„j^+ is hard to extract for n > 12. As 
there are more ways to construct a n-7i^ system, the recursion relation forces us to calculate all 
Qni7c+ ,n27t+ for ^11 pairs ni,n2 = 1---12, before getting to Q1271+.1271+' which costs O(IOO) times 
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Figure 2: The left panel shows energies of a rest system of n-7t^{E„„+) extracted from both methods and 
the right panel shows C1271+ (t), which is decomposed into contributions from excited state, the ground state 
and all thermal states. 
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Figure 3: Left panel shows E„^.^+ „^„;+ extracted from C„^j^+ „^;^+ plotted against x = 2ni +n2- Statistical 
error and systematic error are added up in quadrature. Right panel is C4Qjc+{t) calculated from the new 
method. 



more than the one source case. Similarly studying system of 36 71+ 's requires a third source, 
and becomes 0(100) times more expensive again. In order to overcome this difficulty, we have 
developed a new method which is much faster than the recursion relation method and can easily 
study systems of at least 72 mesons [6]. This new method is based on the fact that the ground state 
energies extracted from all C„,,„2(f) for fixed «i +«2 are the same as they correspond to the same 
spectrum but with different overlaps, shown in Fig. 3. 

For systems containing large number of 71+ 's, thermal states dominate the correlation function 
at later time slices, while excited states dominate earlier time slices, and extracting the ground state 
energy becomes extremely difficulty and systematic error is hard to control. In order to get better 
ground state signals we have also studied the system on 20^ x 256 lattices (with other parameters 
unaltered) using the new method mentioned above. 

3.4 Isospin chemical potential (/i/) and pressure 

The isospin chemical potential, the /X/, is defined as fJ.i{n) = j^, which is approximated by 
using a backward finite difference on the lattice. As the systematic errors of ground state energies 
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Figure 4: The left panel shows energies of a rest system of n-Tt^{E„j^+) extracted from 20^ x 256 lattice, 
and the right panel shows the isospin chemical potential as a functions of isospin density. Lattice space 
a = 0.125 fm. The solid black line is from jPT [2] 



extracted from the T = 128 lattices are large for systems with large number of 71+ 's, we do not 
extract jU/ from these lattices. The 20^ x 256 lattices give much better ground state signals, having 
an extremely long plateau in the effective mass plot, and a one exponential fit is enough to get the 
ground state energy. In Fig. 4, we show the ground state energies and isospin chemical potentials 
from the 20^ x 256 lattice calculated with the newly constructed method mentioned above, /i / is 
consistent with predictions from ;^PT [2] at small density and starts to deviate from them for larger 
densities. 

Since calculations have been done on several lattices sizes {16^, 20'', 24^} x 128, the pressure 
(P) can be derived as function of density by utilizing the discrete version of P = ^ for fixed number 
of 71+. One definition is the following: 



{Li-aY -{Lj-ay 
n 



{{Li + Lj)-a/2y 

where Li,Lj = 16,20,24 and a = 0.125/m. Although we can not extract E„j^+ reliably for large n 
on these ensambles, the pressure can still be studied as a function of isospin density from E„j^+ for 
small n. The pressure is plotted as a function of isospin density in Fig. 5. 



4. Conclusion 

We have calculated correlation functions of systems up to 247r+s in momentum space by ap- 
plying the recursion relation of Ref. [1], verified the dispersion relation, and extracted ground state 
energies of n-n^ systems. Because of the contamination from thermal states in later time slices 
and excited states in earlier time slices, it is extremely difficult to see the ground state of systems 
having more than 1271+ 's for 20^ x 128 anisotropic lattices (a, 0.04) at this pion mass. In order 
to reduce contamination from thermal states, we also use 20^ x 256 lattices. Con^elation functions 
from these larger temporal extent lattices give clear ground state signal, and a one exponential fit is 
sufficient enough to reliably extract the ground state energy even for 72 mesons. 
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Figure 5: The pressure is computed from energy differences for fixed n on different lattice sizes. Red points 
are from 16^-20^ lattices, green points are from 16^-24^ lattices, and blue points are from 20^-24^ lattices. 

Although the recursion relation requires much less time than direct contractions, it becomes 
0(100) times more expensive with an additional source. Studying 3 sources system becomes quite 
expensive and studying 4 sources system becomes impractical. In order to overcome this difficulty, 
we have developed a new method. This new method is far more efficient than the recursion relation 
method, allowing us to study systems of up to 72 mesons [6]. 
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